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ABSTRACT 

The influence of the first stars on the formation of second-generation objects at high redshift may 
be determined, in part, by their metal enrichment of surrounding gas. At a critical metallicity, Zcnt, 
primordial gas cools more efficiently by fine-structure lines of [C II] (157.74 /im), [O I] (63.18 /im, 
145.5 Aim), [Si II] (34.8 /.tm), and [Fe II] (25.99 ^im, 35.35 /im) than by H I or H2 emission. This 
cooling may alter the process of fragmentation into smaller units. We study the time-dependent 
cooling of primordial gas enriched by heavy elements from early massive stars, particularly O, Si, 
and Fe. We define Zdit as the point when the total cooling rate by metals plus H2 equals the 
adiabatic compressional heating. We explore two metallicity scenarios: (1) a single metallicity for 
all heavy elements; (2) individual metallicities {Zq^ Zq, Z^u Z-pc) from theoretical supernova yields. 
For dense gas (n > 10^ cm~"^) with metals in relative solar abundances, fragmentation occurs at 
Zcrit ~ 10^^'^ Zq. However, for lower density gas (71 = 1 — 100 cm^'^), particularly in halos enriched 
in Si, O, and Fe, we find Zcrit ~ 0.1-1% Zq. The critical metallicity approaches a minimum value 
at high-n set by efficient LTE cooling, with thermalized level populations of fine-structure states and 
H2 rotational states (J = 2 and J = 3). Primordial clouds of 10^ Mq at 200 K are detectable 
in redshifted fine-structure lines, with far-infrared fluxes between 10~^^ and 10~^^ W m~^. For 
metallicities Zq « 10~^ and molecular fractions /h2 ~ 10""^, the fine-structure emission lines of [O I], 
[Si II], and [Fe II] could be 10^ - 10^ times stronger than the H2 rotational lines at 28.22 /im (J = 2 — 0) 
and 17.03 /im ( J = 3 - 1). 

Subject headings: cosmology: theory - nuclear reactions, nucleosynthesis, abundances - stars: Popu- 
lation III - galaxies: formation - intergalactic medium 



1. INTRODUCTION 

The last few years have wit nessed an increase in the 
number of studies of first stars llTumlinson fc Shul]|l2000t 
iBromm fc LarsoDl2004l ICiardi k Ferraral2005|) ~7mplicit 

in these studies is the issue of determining the transi- 
tion between the first (metal-free) stars and the second 
(metal-enriched) generation of stars. This transition is 
thought to be determined primarily by the amount of 
heavy elements expelled by early stars and supernovae 
(SNe) and by their rate of incorporation into the sur- 
rounding intergalactic medium (IGM) and gaseous halos. 
At temperatures T > 8000 K, primordial gas cools effi- 
ciently by the excitation of resonance lines of hydrogen 
and helium. To cool be low the "Lyg barrier" at 8000 K, 
the gas must form H2 JPeebles fc Dickelll968t iHirasawal 
^69; LcDD & Shull 198^, whose rotational excitation al- 
lows the gas to cool to lower temperat ures (~ 10^ K) and 
form the first stars ijAbel et alJll997|) . 

First-generation massive stars end their lives as SNe 
and eject newly synthesized heavy elements into inter- 
stellar space and out into the IGM. Once these metals 
are incorporated into neighboring clumps and halos, a 
second generation of stars will form from metal-enriched 
gas. When this gas includes a sufhcient fraction of met- 



als, it can cool more efficiently via the excitation of fine- 
structure levels of key heavy elements (C, O, Si, Fe). 
At this poin t, the gas is said to have reached "critic al 
metallicity" l|Bromm et al]l200lt iBromm fc Loebl l2"0031. 
and the cooling gas is then able to fragment into smaller 
cores and lower-mass stars. 

Metal enrichment and the subsequent cooling both de- 
pend on the mass range of the first stars. Zero-metallicity 
stars above 140 Mq are predicted to end their lives in 
pair- instability supernovae or PISN l|Heger fc WooslevI 
I2OO2I) and eject large amounts of O, Si, and Fe. Massive 
stars in the range 10-100 Mq produce smaller amounts 
of Fe, but synthesize considerable mass in a-process 
elements (O, Ne, Si, S, Ar). As noted by Tumlin- 
son, Venkatesan, fc Shull (2004) and Qian fc Wasser- 
burg (2005), the nucleosynthetic signatures of "Popula- 
tion HI" stars provide key diagnostics of the initial mass 
function (IMF) of the first stars and SNe. 

The purpose of this paper is to calculate a minimum 
metallicity, Zcrit, at which the gas reaches the mass 
scale of fragmen tation after the death of the first gen- 
eration of stars llSchneider et alJl2002l: iBromm fc Loebl 
l2003t iFang fc Cenll2004 . Previous work (|Omukaill2000(l 
used semi-analytical models to describe the thermal and 
chemical evolution of a protostellar cloud with a range of 
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metallicities. In that study, a large network of chemical 
reactions was followed to describe the non-equilibrium 
chemistry for four elements (H, He, C, O) and 45 species, 
but no conclusions were drawn about the critical metal- 
li city. 

iBromm fc Loebl l|2003j) made an analytic calculation 
of the critical metallicity, using cooling rates from fine- 
structure lines of [C II] and [O I]. However, they did not 
include other likely coolants such as [Si II] or [Fe II], 
nor did they follow the thermal evolution of the system. 
Based on the low [Fe/ H] seen in several EMP (extremely 
metal poor) halo stars ijChristlieb et alJ2002() . they iden- 
tified C and O as the elements responsible for the tran- 
sition between Pop HI and Pop II and derived critical 
values, [C/H] = -3.5 and [0/H] = -3.05, for the IGM. 
Their study began with pre-existing clumps of cold, high- 
density gas. They assumed an initial hydrogen density 
uh = 10"* cm~^ and explored a range of temperatures, 
1 00 K < r < 200 K set by Ha cooling. 

iFang fc CenI { 2004) also explored this transition, but 
they used different IMFs and were concerned with the 
photon production of the first objects. To connect Zcrit 
to IGM reionization, they calculated a relationship be- 
tween metallicity of the IGM and ionizing photons from 
PISN and type-II supernovae (SNell). In order to reach 
a critical metallicity of ^ 10~^'^ Zq for an IMF contain- 
ing very massive stars (VMS), the transition between the 
two populations occurred when ~ 1 photon per H atom 
is produced. Approximately 10 photons per H are pro- 
duced with a normal top-heavy IMF. 

In this paper, we assume that a first generation of mas- 
sive stars formed in a high-redshift halo. These stars pro- 
duced a UV background and distributed the products of 
nucleosynthesis throughout the halo. We assume that all 
elements with first ionization potentials less than 13.6 eV 
are kept singly ionized (C II, Si II, Fe II, etc). A simila r 
line of reasoning was followed in iBromm fc Loebl l)2003|) , 
but we introduce three improvements in the astrophysi- 
cal calculations. First, we solve the time-dependent ther- 
mal history of the gas, cooling isobarically from an initial 
temperature and density. Second, we increase the num- 
ber of cooling metal lines, both [C II] and [O I] but also 
[Si II] and [Fe II]. Third, we explore deviations from a 
single metallicity, Z, by defining individual values (Zq, 
Zq, Zgi, Z-pe) in accordance with theoretical yields from 
PISN and VMS. 

The paper is structured as follows. We describe our 
semi-analytical simulations in § (21 We explain our chem- 
ical network and density evolution as a function of the 
different cooling functions in § 12.11 In § 12.31 we give de- 
tails of the metal cooling lines used in the calculations. 
In §121 we present our results. We discuss and summarize 
our conclusions in § 0] 

2. SEMI-ANALYTIC SIMULATIONS 

We have approached the cooling problem by consid- 
ering a parcel of gas with primordial composition (see 



Appendix). We give the gas an initial temperature 
T w 100 — 200 K and initial metallicity and let it cool 
isobarically, until the gas reaches a condition for frag- 
mentation: a rate of cooling faster than adiabatic heat- 
ing rate within a local Hubble time. We adopt an initial 
mass density of hydrogen, pvir ~ 178/3b(z), as that of 
a virialized halo of total mass M — (10^ Mq)Mq and 
baryon mass {yib/Q,m)M at fiducial redshift z k, 20. 
The hydrogen number density in all forms is given by 
nn = (1 — Y)p/mH, for helium abundance Y — 0.244 by 
mass and /3b^o=4.21xl0^^^g cm~^, appropriate for the 
baryon dens ity parameter Uhh"^ = 0.0224 ± 0.0009 found 
by WMAP (jSoergel et al.ll2003|) . The virial radius and 
average hydrogen density of the halo are then, 
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Gas within this halo probably contains clumps of 
higher density, but this formula provides a characteris- 
tic mean value. This mean density is 178 times larger 
than the mean density of the IGM, (nn) = (1.90 x 



10- 



70 



We chose a value of IStt^ « 178 



for the density ratio at virialization, its value for il = 1. 
For a cosmology with fl + A = 1 and = 0.3, this ratio 
is ~ 340. In general, the total matter density is 



(3) 



where Ac(ri, A) w ISTr^f^" " ^^, a formula ac curate to 5% 
for values of 17 = 0.15-1.0 CEke et al.lll998|l . 

The condition for fragmentation has been investigated 
by many authors, both analytically and through hy- 
drodynamical numerica l simulations (Rees fc Ostrikeij 
1977t iBlumenthal et ab 1984: Yama da fc Nishi 1998 
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tion, tcooi tg, requires that the cooling time be shorter 
than the free-fall time, the time for a perturbation at 
the surface to reach the center of the clump, 
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Here, rit is the total particle number density and C{n, T) 
is the cooling rate per unit volume (erg cm~'^ s~^). In 
this cooling rate, we include [C II], [Si II], [O I] and [Fe II] 
fine-structure and forbidden lines, as well as cooling by 
H I and H2. These four heavy elements (C, O, Si, Fe) 
are products of a-burning whose dominant interstellar 
ions have efficient cooling from low-lying fine-structure 
levels. Other a-process species such as Ne I, Mg II, S II, 
Ar I, and Ca II lack fine structure states in their ground 
configurations. We refer the reader to § 12.21 and § 12.31 
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for a complete description of the atomic data and metal 
cooling rates. 

As the cloud contracts, with sufficient cooling for fra- 
mentation, the temperature drops and the density in- 
creases. The cloud sizes at fragmentation are character- 
ized by the Jeans length scale and Jeans mass, Mj. Frag- 
mentation slows when the gas becomes optically thick, 
but also when the gas density exceeds the critical densi- 
ties for coUisional de-excitation of the dominant coolants. 
This defines a minimum mass, which we use as the con- 
dition for fragmentation. This c ondition depends on the 
choice of the cooling function. iBromm fc Loebl l|2003f) 
studied the formation of low-mass stars from gas that has 
been "polluted" by trace amounts of carbon and oxygen. 
To derive the critical metal abundances for fragmenta- 
tion, they looked for the point at which the metal- only 
cooling rate, >Ccii,oi, equals the adiabatic compressional 
heating rate, Fad- 

The gas fragments when C = Fad, where Fad = 
—P{d \np/dt) = 2>P{dR/dt)/R is the adiabatic compres- 
sion heating per unit volume (erg cm~'^ arising from 
free-fall collapse of the gas under gravity. 

Fad = (247r)i/2 [1 _ R/Ro]'/\n,kT)^ . (5) 

Here, p = 1.323n//m// is the total mass density, where 
uh is the hydrogen number density and rit = LOSOlnn is 
the total particle number density (mostly H and He, since 
the electron and molecular fractions are quite small). As 
before, these values assume helium mass fraction Y = 
0.244 (He/H = 0.0807 by number). 

The initial and current cloud radii are Rq and i?, and 
the coefficient (247r)i/2 [1 - i?/i?o]^/^ varies from 1.22 to 
2.74 for clouds with radius (R/Rq) =0.99-0.90. We use 
the value corresponding to a 10% change in radius, for a 
maximum coefficient (2.74), and we write the adiabatic 
heating rate (erg cm^'^ s~^) as. 



Fad = (4.59 X 10-^") 
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The different terms that form part of the cooling func- 
tion used in these definitions are explained in the follow- 
ing section. We use the fragmentation criterion C — Fad 
and explore cases in which adiabatic heating is compared 
to metal-line cooling only, as well as to the total cooling. 
In practice, both methods give similar results. 

2.1. Chemistry network and thermal evolution 

Throughout this paper we assume a ACDM cosmology 
with fio = 0.30, riA = 0.70, f^bo = 0.0457, h = 0.7, 
and CTg = 0.90, the root-mean-square density disper- 
sion within spheres of radius 8/i^^Mpc. As initial con- 
ditions, we adopt residual electron and molecular frac- 
tions ricrcs/'T-H = 5 X 10~^ jScag er et al. 2000) and 
riH2.ros/nH = 2 X 10"-* ()Lepp et al.ll2002() . These frac- 
tions arc "frozen out" from expansion following the re- 
combination epoch. We model the gas as 75.6% H and 
24.4% He by mass and include traces amount of metals 



(7) 



as coolants. We follow the non-equilibrium chemistry of 
nine species, Hz, H, H"* ", HJ, H^, He, He+, He+^, and e~ 
()Hutchings et alJl2002(l . The reaction rates and chemical 
processes can be found in the Appendix. 

We follow the evolution of the gas density at constant 
pressure, P. From the energy equation, we find, 

dnt nt C 

~dt " P~ ' 
where rit is the total particle number density and C is 
the cooling rate per volume. From the ideal gas relation, 
we obtain the temperature, T = P/utk, and define 

^ = ^metals + -Chs + /^H + 'Cho+ , (8) 

where the individual terms denote the cooling rates from 
metals, H2, H°, and He+. 

To find the molecular hydrogen cooling function, we 
solve for the statistical equilibrium of the first six rota- 
tional levels (J = 0-5). We use the radiative transition 
rates for the H2 rotational quadrupole lines from Wol- 
niewicz, Simbotin, & Dalgarno (1998): A{J = 2^0) = 
2.94 X 10-" s-i, A{J = 3^1)= 4.76 x s^^, 
A{J = 4^2) = 2.75 X 10-9 s^^, and A{J = 5 
3) = 9.83 X 10^9 s^^. We assume separate populations 
of ortho-Ha (J = 1, 3, 5) and para-Ha (J = 0, 2, 4) with 
the usual 3:1 ortho/para ratio from nuclear spin multi- 
plicities. The rate coefficients for H° and H e° collisional 
excita tion of these levels were taken from iForrev et alJ 
1I1997I) . who calculated new rates of (non-reactive) rota- 
tional excitation over the range 10 K < T < 1000 K, 
using exact vibrational wavefun ctions and a more reli- 
able H3 potential energy surface ( Boothrovd et aljri996() 
for the interaction potential . Their H°-H2 rates ar e sev- 
eral times larger than those IjLe Bourlot et al .I1999I) com- 



monly used in other studies. We include He-H2 col- 
lisions, whose rate coefficients are 30-40 times larger 
than those of H°. We scale the He rates to H rates of 
IForrev et alJ ( 19971). using ratio s of He/H rate coefficients 
computed by Le Bourlot et al.1 (^999) . 

For the complete set of collisional excitation rates 
and other less important p rocesses together with their 
analytical fits see IHutching s et al ] §002). The atomic 
data for the metal fine-structure and forbidden lines 
and the cooling rates, /Imctais, are explained in the 
following sub-sections. Table ^ shows the r eference 
solar m etallicities used in this work, taken from iMortonI 
()2003|) for C, O, Si, and Fe. The table also fists 
the critical number densities, ricr = ^21/721, for H° 
de-excitation of the fine-structure lines, where the level 
populations approach Boltzmann (LTE) populations. 
Here, A21 and 721 are the radiative decay rate and 
H°-impact de-excitation rate coefficient respectively. 
For n > Tier, the volume rate of line cooling scales as 
C oc nZcrit, rather than n^Zcrit at low density. Because 
Fad OC Tn^/'^, the critical metalHcity Zcrit 
exhibits an approximately parabolic shape around its 
minimum at n w ricr for a particular coolant. For a 
given heavy element, Z^rsx scales as n~^l'^ for n < ncr. 
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TABLE 1 

Metals and Transitions Used in Calculations^ 



Species 


A(/im) 


{Ai/AH)Q 


(Ai/Ape)© 


ncr(cm 


C II 


157.74 


2.45(-4) 


7.75 


2.86(3) 


O I 


63.18 


4.90(-4) 


15.5 


6.08(5) 


Si II 


34.8 


3.63(-5) 


1.15 


2.75(5) 


Fe II 


25.99,35.35 


3.16(~5) 


1.00 


2.24(6), 1.46(6) 



''Notation: 2.45(-4) denotes 2.45 X 10"*, A is the wavelength of 
the fine-structure transition, and {Ai/ Ah)i3 is the solar abundance 
of element (i) relative to hydrogen (Morton 2003). 

^ricr is the critical density of the transition, at which H°-impact 
collisional de-excitation (T = 200 K) equals radiative decay. 



and as n^/^ for n > ricr- The most efficient cooling 
occurs at n « ricr for the fine-structure hnes of a given 
element. 

2.2. Atomic Data for Metal Lines 

2.2.1. C II fine- structure line 

The ground state (2p) configuration of C 11 is treated 
as a two-level atom, with a fine-structure radiative tran- 
sition from (2p) -^"^ Pm-)]- The atomic parameters 
were taken from lHollenbach fc McKeel l|f 989l) . This tran- 
sition corresponds to 157.74 /zm with excitation energy 
E21 = 1.259 X 10"^'* erg and excitation temperature, 
Tcxc — E2i/k — 91.2 K. The radiative transition rate is 
A21 — 2.4 X 10^^ s^^, and the H°-impact de-excitation 



rate coefficient is: 

72^1 = (8.0 X 10-1" cm^ s-')T?oo , (9) 

where Tioo = (T/lOO K). The critical density for H° de- 
excitation is ncr = A21/721 = (3000 cm-3^'T^-"•07 
cooling rate is 



J . 1 ue 



CII 



(10) 



72" +71^2 + ^21 /nn ' 

where Zq is the carbon metallicity and Ac is the so- 
lar carbon abundance relative to hydrogen. By detailed 
balance, we have, 

7" = (.92/31)72" exp{-E2i/kT) , (11) 
where 52 = 4 and gi — 2 are statistical weights. 

2.2.2. O I fine- structure and forbidden lines 

The ground-state (2p^) configuration of O I is a 5- 
level multiplet, with transitions between fine structure 
levels arising from (3P2,3Pi,3Po,iD2, ^Sq). We label 
these levels, from lowest to highest, as 1-5; their ex- 
citation temperatures from the ground ('^P2) state are 
Texc = 227.7 K (3Pi) and 326.7 K (^Po). For the ex- 
cited terms, Texc = 22,523 K (ID2) and 48,661 K (^Sq). 
The fine-structure transitions within the ^P term lie at 



were taken from lHollenbach fc McKe'3 



far-infrared wavelenegths: 63.18 fj,m (J — 1 —^ 2) and 
145.53 fiui (J = — > 1). The optical forbidden lines lie 
at 6391.5 A, 6363.8 A, and 6300.3 A (from to ^P) and 
at 2972.3 A, 2958.4 A, and 5577.4 A from the term. 

The radiati ve transition proba bilities were taken from 
tabulations in lOsterbroc^ ljl989|) : A21 = 8.9 x 10"^ s^^, 
A31 = 1.0 X 10-1° s-\ A32 = 1.7 X 10-5 s-\ Ail ^ 
6.3 X 10-^ s-i, A42 = 2.1 X 10-3 s-\ ^43 = 7.3 x 
10-^ s-i, A51 = 2.9 X 10-1 s-i, A52 = 7.3 X 10-^ s-\ 
^54 = 1.2 s-i. The H° collisional de-excitation rate coef- 
ficients (cm'^ s~ 
l|l989(l : 7^ 

7II = 1.1 X 10-i°T°o^''. The remainder of the rate co- 
efficients are less well determined and are estimated to 
be 10-1^ cm^ . The critical density of the first ex- 
cited state, 3Pi by H° de-excitation is ncr = (9.67 x 
10^ cm-3)rK,o^'^. The radiative cooling rate from [O I] 
was computed from the equilibrium populations, multi- 
plied by their radiative decay rates and line energies. 

2.2.3. Si II fine- structure line 

For the Si II (3p) ground-state configuration, we also 
provide a two-level treatment, with a fine-structure ra- 



9.2 X 10-iirOoe^ 73^ 



4.3 X lO-iirfoSO, 



diative transition from (3p) [ 



3/2 



Pi 



/2J 



This tran- 



sition corresponds to 34.8 //m with E21 — 5.71 x 10-^^ 
erg and Tb x^ = 410 K. The atomic p arameters were 
taken from iHoUenbach fc McKeel l)l989() . The sponta- 
neous transition rate is A21 = 2.1 x 10-^ s-^, and the 
II°-impact de-excitation rate coefficient is: 

-0.07 



72^1 



(8.0 X 10-^" cm^ )Two ■ (12) 
The critical density for 11° de-excitation is = (2.62 x 
)T°5q^. As in equation 1101 the cooling rate is: 

7}^2^2i-E2i ^'si^si n-a 
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Sill 



2.2.4. 



721 



-7?2 



-42l/nH 

Fe II fine- structure lines 



(13) 



The ground-state (^D) configuration of Fe II is mod- 
eled as a 5-level atom, with fine-structure levels of to- 
tal angular momentum J = 9/2, 7/2, 5/2, 3/2, and 1/2 
(labeled 1-5 respectively). The fine-st ructure level ener- 
gies were taken from the NBS tables ()Moorelll952|) . and 
the radiative transition rat es were ado pted from calcula- 
tions of the "Iron Project" l)Quinet et al. 1996). The fine- 
structure levels 2-5 have excitation temperatures Texc — 
553.6 K, 960.6 K, 1241.1 K, and 1405.7 K above ground. 
The primary radiative transitions occur at mid-infrared 
wavelengths: 25.99 /im (J = 7/2 9/2), 35.35 /.tm 
(J = 5/2 7/2), 51.28 ^m (J = 3/2 5/2), and 
87.41 //m (J — 1/2 3/2). The spontaneous transition 



rates are: A2 



2.13x10-'^ s" 



A 



32 



1.57 x 10-3 s- 



X 10" 



taken 



1.50 X 10-9 s-\ 



^43 = 7.18 X 10-^ s-i, = 1. 
from lOuinet et"al] lll996l), and A :> 
taken from iHollenba^^lTfc McKeel I 

The H°-impact de-excitation rate coefficients were 
taken from .HoUenbach fc M cKcc (1989), with no as- 
sumed temperature dependence. These rates (in units 
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are: 7^ = 9.5 x IQ-^^ 7^ = 5.7 x 10-^° 
and — 4.7 x 10^^". In the absence of informa- 
tion on higher-J levels in the multiplet, we have 
adopted de-excitation rates of 5 x 10^^" cm'^ s~^, con- 
stant with temperature. The critical density of the first 
excited state, ^07/2, resulting from H° de-excitation is 
ricr = 2.24 x 10^ cm~^. We computed the radiative cool- 
ing rate from [Fe II] from the statistical equilibrium pop- 
ulations of the excited fine-structure levels, multiplied by 
their radiative decay rates and line energies. 

2.3. Metal cooling 

At T < 8000 K, excitation of Lya is no longer an ef- 
ficient cooling process. Similarly, at T < 60 K, cool- 
ing by pure rotational lines from H2 becomes ineffec- 
tive. At very low temperatures (T < 60 K, HD cool- 
ing can be effective, despite the small primordial abun- 
dance, (D/ H)nrini ~ 26 ppm, owing to its weak d ipole 
transitions l|Lepp k Shul]| ll984t lEepp et alJl2002(l . At 
T < 1000 K, in the absence of external photoionizing 
sources, the hydrogen gas is completly neutral, and low- 
energy fine-structure transitions of [C II], [O I], [Si II], 
and [Fe II], excited by H° collisions, become the most im- 
portant coolants. (We assume that C, Si, and Fe are pho- 
toionized by ambient UV radiation from newly formed 
hot stars.) These fine-structure transitions can also be 
excited by electrons, once the electron fraction, Xe-, rises 
above about 10"'^. At redshifts, z < 20, the residual 
electron fraction is Xe ~ lO"'^'^, and the electrons con- 
tributed by photoionization of trace metals (0+, Si"*", S"*", 
Fe+) are negligible. As the cloud cools and condenses, 
Xe declines owing to recombination of the trace H+ . 

In Figures ^ and El we plot the four metal cooling 
functions together with £h2 at metallicities of 10~^Zq 
and 1O~^Z0, respectively. In both figures, we show the 
total cooling function, as well as individual rates from 
H2, [C II], [O I], [Si II], and [Fe II]. The low-density 
cooling rates are for gas with total number density 
nt = 10~^ cm~^. At high densities 
and low metallicities {Z < 10~^ Z, 
excitation of H2 becomes important at T < 200 
compared to the metal lines. 



nH > 10 cm 
o), cooling by H° 
-400 K 



2.4. Analytic Estimates of Critical Metallicities 

The criterion for fragmentation, used to derive Zcnt, 
comes from equating heating and cooling functions of 
the form. Fad = aTrt'/'^ and C{n,T). From equation (5), 
we have a = 2.74fc(G/i)^/^, where /i — 1.323m//. Above 
the critical density, when the fine-structure levels reach 
LTE (Boltzmann) populations, the LTE cooling function 
can be written £ — nZf{T), where /(T) can be derived 
from the LTE excited populations times the radiative 
transition rates (^21) and line energies (i?2i)- Thus, the 
critical metallicity is 
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(14) 



c in 

<ll 

\ I 



-1 — I — 1 — I- 



-l 1 — 1 — 1 — 1 — I — I- 



ni=10"'' cm 




I 



1 2 
logio(T/K) 



cn 
o 



1 1 1 
■ ni=1000 cm"' 


1 1 i 1 1 1 1 1 


i i 1 i 


- z=io-% 










■-■'H2_ _ - 




^ . ' ^ ' — ^ 


- ' " " Fell " 


cm/ 


/ /• / y 

' : :/ 
! :' '1 
1 ■■ .■/ 
0]- .sill 

i // 






/ .■ ,7 






11'/ 


1 



I 



1 2 
log,o(T/K) 



Fig. 1.— Total cooling rate at Z = 0.01 Zq (thin full line) 
and individual rates for [C II] (thick full line), H2 (dash- three- 
dotted line), [O I] (dash-dotted line), [Si II] (dotted line) and [Fe II] 
(dashed line). The gas is assumed to have total number densities 
rit = 10~^ cm~^ (top panel) and 10^ cm~'^ (bottom panel). 



The LTE cooling function, /(T), can easily be evaluated 
for a two-level atom of element (i) of solar abundance Ai 
relative to hydrogen: 



Jii ^21 -^21 



[l + (gi/52)exp(4lVfcT)] 



(15) 



where gi and 32 are the statistical weights of the lower 
and upper states, respectively, of the two-level transition. 
From the above, we arrive at a convenient analytic ex- 
pression for the critical metallicity for a given atom and 
fine-structure coolant, that depends only on its reference 
abundance and the atomic parameters of the transition: 
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1 -I- — exp{E^^/kT) n^/^T 
(16) 

We can evaluate this expression for the two-level species 
[C II] and [Si II] at n > n^,- and T = (200 K)T2oo, where 
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Fig. 2. — Cooling rates, as in Fig. 1, for gas with metallicity 
Z = 10~*Zq and densities nt = 10~^ cm~^ (top panel) and 10'^ 
cm~'^ (bottom panel). Metals have abundances in solar ratios. 



efficient fine-structure LTE cooling begins: 



4^,"^ = (3.74xlO-'^Zo)(^ 



n 

ricv 



,(SiII) 



.34 X IQ-^Zq) 



1/2 



1/2 



(17) 



(18) 



For [O I] and [Fe II], with 3-level and 5-level fine- 
structure states, respectively, we can apply the 2-level 
formalism to the lowest two states, to find approximate 
values, 



7(01) ~ / 



^(1.02x10-^Zq)( — 

flcr 



zir^«(1.76xlO-4Zo)(^ 



1/2 



1/2 



(19) 



(20) 



The critical densities for H° de-excitation are: ncr(CII) = 
(2860 cm-3)r2^°°^ n„(OI) = (6.08 x 10^ aji-^)T-°Q^\ 
ncr(Sill) = (2.75 x 10^ cui-^)T°q°'^ , and ncr(Fell) = 
(2.24 x 10^ cm^'^). These formulae provide reference 
metallicities, but they are no substitute for the exact 



calculations described below. For instance, they do not 
give the temperature dependence implicit in equation 
(16), arising from the Boltzmann factor for excitation 
of fine-structure states. However, for LTE conditions at 
high density, these formulae agree fairly well with the 
numerical computations. Small differences occur for the 
multi-level [Fe II] and [O I] lines, since we used a 2-level 
approximation that becomes inaccurate at higher tem- 
peratures by omitting higher states. 

3. RESULTS 

In S I3.1l we show results from simulations with a single 
metallicity, Z, adopting the same relative metal abun- 
dances for each of the four metals considered. In S l3.2l we 
include individual metallicities for C, O, Si, and Fc, using 
metal yields from PISN calculated by iTumlinson et alJ 
()2004|) . In our formulation, each metal appears with a 
different abundance relative to the other three elements, 
depending on the mass range of progenitors. Our anal- 
ysis of Zcrit follows the fragmentation criterion (§ 2), 
in which we compare the adiabatic heating with cool- 
ing rates from metals plus II2. We also present the final 
properties of the fragmented gas, such as temperature 
and Jeans mass, for a variety of metallicities and densi- 
ties. 

3.1. Single-metallicity simulations 

First, we assume that all metals in the gas appear with 
the same relative abundance, characterized by a single 
metallicity (Z). Figure synthesizes the results of many 
models of cloud evolution, for different initial densities, 
with an initial temperature of 200 K. The elements are 
color-coded in groups (C-alone, O plus Si, Fe-alone) re- 
flecting their probable origin from progenitor mass ranges 
and nucleosynthetic pathways. There are three set of 
curves corresponding to three choices of the cooling func- 
tion. In red, we use only metal cooling by [C II] . In green, 
we use the combined metal cooling function of [Si II] and 
[O I], and black corresponds to only [Fe II]. We have 
combined the results for [Si II] and [O I] , since these ele- 
ments come from similar massive stars, and because this 
presentation makes the results less cluttered. The [O I] 
63 fim line dominates the cooling when only those two 
elements are considered. 

Each set of three curves (red, green, black) in Fig- 
ure 01 corresponds to the locus of fragmentation points 
in metallicity Z and final density n/. The code starts 
with an initial density, corresponding to the point far- 
thest to the left of each curve. Each point on the 10 sets 
of curves is found by running many simulations for var- 
ious initial densities and metallicities, until they satisfy 
the fragmentation criterion (C = Fad) at Zcrit and den- 
sity n / . The double- valuedness in many of the parabolic 
curves is explained by Figure 0] which shows the cooling 
rates and equilibria for gas with Ui ~ \ cm^"^ and three 
metallicities. 

Table 121 summarizes the results of Figure by choosing 



7 




TABLE 2 
Critical Metallicities'^ (log Z^it) 



Fig. 3. — Critical metallicity as a function of final number den- 
sity. Each parabola shows the density evolution for a gas for 10 
different initial number densities. From left to right: 10~^ cm"'' 
to lO'^ cm~'^ in order-of-magnitude steps. Red corresponds to gas 
composed only of hydrogen plus trace C II, green shows Si II and 
O I, and black shows Fe II. 




1.5 
log,o(T/K) 

Fig. 4. — Radiative cooling rates and equilibria for gas with 
rii = 1 cm~'^ and three metallicities. The adiabatic heating rate is 
also shown. Note the double solutions at log T = 1.5 and 2.2 for 
Z = 0.007 Zq. These two thermal equilibria are reflected in the 
parabolas in Figure 3. There are no equilibria at Z = 0.003 Zq, 
so that fragmentation does not occur. 



three particular values of temperature and density rele- 
vant for this study. The first column shows the tempera- 
ture and density of the gas. For the low-density case, we 
choose the set of three curves corresponding to an initial 
density of 0.1 cm~^, while the high-density case starts at 
10** cm~^. We begin at = 200 K, and after the gas 
cools down, we choose two more temperatures, 150 K 
and 100 K and their corresponding densities. The next 
four columns show the metallicity at those temperatures 
and densities during fragmentation. In the fourth col- 
umn, all metal lines are included in the cooling function. 
The last column shows the Jeans mass, 



(T, n) 


C 11*= 


(O I, Si II) 


Fe II 


All Metals 




Low n 


200, 1.0{-1) 
150, 1.3{-1) 
100, 2.0{-l) 


-1.59 
-1.69 

-1.80 


-1.01 
-0.92 
-0.67 


-0.16 


-1.70 
-1.77 
-1.83 


2.1(5) 
1.2(5) 
5.7(4) 


High n 


200, 1.0(4) 
150, 1.3(4) 
100, 2.0(4) 


-3.22 
-3.20 
-3.19 


-3.58 
-3.44 
-3.16 


-2.75 
-2.46 
-1.86 


-3.78 
-3.66 
-3.49 


698 
353 
135 



''Minimum values of log Zcrit together with final values of density 
n (in cm"'') and temperature T (K) at the point of fragmentation. 
For [Fe II], we found no sensible low-Z solutions at T = 100 K and 
150 K, since its first cooling level lies at Texc = 554 K. 

''Columns 2-5 give the critical metallicity, in units of solar metal- 
licity, corresponding to a choice of metals in the cooling function. 
Column C II means that only [C II] lines are included; Column 
labeled "All Metals" includes C, O, Si, Fe. 

Jeans mass Mj in Mq. Same format as in TablelTI where 1.0(4) 
denotes 1.0 X 10"*. 



at fragmentation, 
Mj = (700 Mq) 



3/2 



( ^ 

V104 cm-3/ 



-1/2 



(21) 



200 Ky Vl04cm-3- 

taken from lClarke fc Bromml l)2003j) . where T^ag is the 
temperature at fragmentation and hh is the hydrogen 
number density. 

The minimum critical metallicities are calculated from 
Figure |31 and shown in Figure |31 The bottom panel of 
FigurelHshows two sets of curves, to illustrate the sHght 
differences in fragmentation criteria when one equates 
Fad to cooling by metals-only (dashed lines) and to to- 
tal cooling (solid lines). In both cases, one can see the 
transition from non-LTE to LTE of particular emission 
lines. The minimum value in each curve corresponds to 
the critical density of the line (see Table QJ. The curves 
turn around when coUisional de-excitation of the lines 
starts to dominate for each metal. The volume cooling 
rate then scales as n rather than n^. The cooling effi- 
ciency is reduced, and higher abundances of metals are 
needed to reach fragmentation. 

Table|3|shows the values of the minimum critical metal- 
licities for the fragmentation criterion, and for several 
values of gas density. We choose a range in densities, 
from n = 0.1 cm"'^, just below the mean density of gas 
in virialized halos at z = 20, up to n = 10"* cm~^, charac- 
teristic of dense cores. Within these halos, and particu- 
larly at their center, one expects to find clumps of cooler 
gas near the threshold of gravitational instability. These 
star-forming cores may have cooled or been compressed 
to higher values. There may also be multiple clumps 
within the halo. These inhomogeneities probably play 




TABLE 3 
Minimum Critical Metallicities^ 
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Fig. 5. — Top panel. Minimum critical metallicities for individual 
heavy elements vs. total number density at T = 200 K. Curves 
correspond to gas enriched only by C II (dashed), Si II (dash-dot), 
O I (dash-triple-dot), and Fe II (dotted). Solid curve shows all 
four species together. Bottom panel. Solid lines are simulations 
that use a fragmentation criterion that equates adiabatic heating 
to cooling by metals only, while dotted lines assume cooling from 
metals plus H2 (see §0. Minima occur at high densities, with 
log Z^rit = -3.48 (C II), -3.54 (Si II), -3.78 (O I), -3.52 (Fc II), 
and —4.08 (all elements). 



an important role in how the first star formation is initi- 
ated. They may also affect how the heavy elements are 
incorporated into surrounding gas, how those gas clouds 
are taken to critical metallicity, and how subsequent gen- 
erations of stars are produced. 

Figure shows the final state of the gas at the in- 
stant of fragmentation. We fix the metallicity and evolve 
the cooling of the gas, for many initial densities, from 
T — 200 K until the fragmentation criterion is satis- 
fied. The final temperature and Jeans mass of the frag- 
ments are plotted as a function of the final density of 
the gas. In both panels, we ran the simulation for five 
different metallicities, including all the metal lines. The 
curves correspond to metallicities ranging from 10"'* Zq 
to 1 Zq. Larger amounts of metals allow the gas to 



rii 


C II 


(O I, Si II) 


Fe II 


All Metals 


(cm~^) lo 


g -^crlt 


log Zcrit 


Og .Zcrit 


log Zcrit 


0.1 


-1.77 


-1.00 


-0.16 


-1.86 


1.0 


-2.35 


-1.53 


-0.69 


-2.36 


1.0(1) - 


-2.84 


-2.07 


-1.23 


-2.87 


1.0(2) - 


-3.29 


-2.57 


-1.73 


-3.31 


1.0(3) - 


-3.38 


-3.15 


-2.29 


-3.67 


1.0(4) - 


-3.23 


-3.59 


-2.75 


-3.79 


1.0(5) - 


-2.72 


-3.91 


-3.15 


-4.01 


1.0(6) - 


-2.22 


-3.92 


-3.49 


-4.05 



''Minimum values of log Z^rit needed for fragmentation of gas with 
initial temperature T = 200 K and number densities rii . 



cool to very low temperatures, as expected. Figure |B1 
also illustrates the possibility of thermal phase transi- 
tions, resulting in a drop in the final temperature. The 
temperature reaches a minimum at rij^ « lO'^'^ cm~'^ for 
-3.3 < log [Z/Zq) < 0, but rises to n/ « 10^ cm'^ 
for Z < lO"^'*^ Zq. This shift results from a change in 
primary coolants, from [C II] cooling {ricr ~ 3000 cm^"^) 
to fine-structure lines ([O I], [Si II], [Fe II]) with higher 
critical densities, 3 x 10^ cm"-^ to 2 x 10^ cm"'^. 

Figure [T] shows the same final temperature as in the 
top panel of Figure |H| emphasizing the effects of specific 
heavy elements. As labeled, wc show cases with no metal 
lines, only Fe II, only Si II plus O I, only C II, and all 
four elements in the cooling function. In each case, the 
metallicity is Z = 3 x IQ-'^Zq. At n w 2 x 10^ cm'^, the 
cooling is driven mainly by [O I] lines. This is because 
the C II levels reach LTE, so its cooling scales as n rather 
than n^. 

3.2. Multiple-Metallicity Simulations 

Empirical evidence from extremely metal poor (EMP) 
halo stars (Christlicb ct al. 200.2^) suggests that many 
progenitor stars, with a range of masses, contributed to 
the observed nucleosynthetic patterns. This point was 
stressed by Tumlinson, Venkatesan, & Shull (2004) and 
Qian & Wasserburg (2005). When exploring the concept 
of "critical metallicity" , we need to understand the ef- 
fects of different mixtures of heavy elements, produced 
by massive stars and SNe with a range of initial masses. 

We adopt the relative abundances of metals from stud- 
ies o f PISN yields. More sp ecifically, we use the results 
from lTumlinson et al.l ( 2004) who reviewed and critiqued 
the hypothesis that the primordial stars were VMS. They 
calculated the theoretical yields from PISN for three pro- 
genitor mass ranges: 150-170 Mq, 185-205 Mq, and 
240-260 Mq (see their Figure 1). We use these yields 
and run three sets of simulations, in which we take the 
Fe abundance as a free parameter. Depending on the 
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Fig. 6. — Top panel: fragmentation temperature as a function of 
the final density, rif, for gas with initial temperature T = 200K. 
Bottom panel: the Jeans mass at the point (Tfjag, fifrag) in units 
of solar masses. All curves are labeled by the metallicity, ranging 
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Fig. 7. — Fragmentation temperature for a gas with no metals 
(red line), only Fe II (dark blue line), only Si II plus O I (green 
line), and only C II (black line). The light blue includes all four 
metals. In all cases, the metallicity is fixed at Z = 3 X 
and the initial gas temperature is T = 200 K. 



TABLE 4 
Metallicity yields from PISN^ 



Mass Range [C/Fe] [O/Fe] [Si/Fe] [Si/O] 



150 - 170 Mq 
185 - 205 Mq 
240 - 260 Mq 



1.4 
-0.2 
-1.3 



1.9 

0.4 
-0.8 



2.4 

1.2 
-2.0 



0.5 
0.8 
-1.2 



Metalli city yields from PISN iHeger fc WooslevI 120021) as inte- 



grated bv lTumlinson et al.l ^20n4^ 



range of masses and their yields, we calculate the metal- 
licity of the other three elements (C, O, Si). Table 0] 
gives the abundances of [C/Fe], [O/Fe], and [Si/Fe] and 
the mass ranges. 

Figure |H| shows the minimum metallicity needed for 
fragmentation for three mass ranges of PISN. The solid 
line in all panels is the critical metallicity. For each mass 
range, Zait varies considerably, depending on the dom- 
inant elements synthesized. For a mass range 150-170 
Mq , the main component of total Zcrit comes from Si II, 
and the less important component is Fc II. For the high- 
est mass range, 240-260 Mq, the main metal cooling 
comes from Fe II, while the least important element is 
Si II. To a good approximation, C II and O I keep the 
same relative importance in all three mass ranges, with 
O I contributing slightly more cooling than C II. 

Figure|n|summarizes all the results. The full line shows 
Zcrit, considering a single metallicity, while the dashed 
and dotted lines show Zcrit from simulations with indi- 
vidual metallicities of C, O, Si, and Fe, given by the 
theoretical yields of PISN. 

3.3. Luminosities and Detection of Cooling Lines 

In the preceding sections, we have explored the mi- 
crophysics of fine-structure cooling of primordial g as, as 
it makes a transition from H2 cooling to metal fine- 
structure line cooling. An obvious question is whether 
such cooling clumps are detectable at high redshift. How 
much mass must be cooling, and to what metallicity must 
it be enriched, in order to rise above the threshold of de- 
tectability by current and future telescopes? 

For redshifts in the range 3 < z < 10, the fine-structure 
lines under consideration will shift into the far-infrared 
(FIR) and sub-millimeter regions of the spectrum. On 
the ground, partial-transmission windows exist in the 
sub-mm, near 350 /xm (333-380 ^m), 450 nm (420-500 
^m), and 850 ^m (830-910 ^m). These FIR windows 
will be observable on the ground with ALMA. In space, 
there are several planned missions that offer a chance 
for more sensitive FIR observations with mirrors cooled 
to 5-10 K. The SPICA mission is a proposed Japa nese 
infrared astronomical satellite I Matsumotcl l2003|) de- 
signed to explore the universe with a cooled (4.5 K) tele- 
scope of 3.5m aperture. A spectrometer aboard SPICA 
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2 4 

Fig. 8. — Minimum critical metallicity for fragmentation. Each 
plot corresponds to one of the three mass ranges of PISN yields. 
The solid lines show the minimum critical metallicity, while the 
dotted and dashed lines are the components of the solid line, [C II], 
[O I], [Si II], [Fe II]. Top panel gives 150-170 M0 mass range, 
middle panel shows 185-205 Mq range, and bottom panel shows 
240-260 Mq range. In the top panel, the dominant cooling is from 
[Si II] , and in the bottom panel it is from [Fe II] . These shifts result 
in changes in .Z^ritC") relations, since Fe-rich gas has higher critical 
metallicity (see Fig. 5). 
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Fig. 9. — This plot summarizes the two metallicity assumptions 
used in this paper. Full line shows the total minimum critical 
metallicity from Figure |F] All four elements appear in the cooling 
function with the same relative abundances. Dashed and dotted 
lines show the same results, but with different relative abundances 
given by the yields of three mass ranges of PISN. 



( Bradford et alJl2003j) could reach 5cr flux limits of lO^^i 
W m~^. An even more ambitious mission is SAFIR (5*171- 
gle Aperture Far Infrared Telescope), a lOm-class, cold 
(4-10 K) space telescope with hopes of achieving broad- 
band FI R sensitivity to emission-line fl uxes of 10^^^ 
W m-2 llLesteill2004l iBenford et alll2004^ . 

To produce detectable FIR line emission, the cooling 
gas clouds must have baryon masses of at least Mb — 
(10* Mq)Ms at metalhcities Z > 0.01 Zq. If the clouds 
have high densities, nn > f^cr, the cooling rate is propor- 
tional to mass, and we can estimate the line luminosities 
from the LTE cooUng rates, Li = c'^^^{T){Mb/ n). Us- 
ing the results of § 2.4 for T w 200 K, we estimate that 
the metal fine structure lines have luminosities: 

Leu = (3.7 X 10^9 erg s-i)M8(Zc/0.01 Zq) (22) 



-^Sill - 



1 X lO''" erg s~i)Af8(^si/0.01 Zq) (23) 



Loi = (2.0 X lO-^i erg s-^)Ms{Zo/0.01 Zq) (24) 
LfcII = (2.2 X 10^1 erg s-1)A/8(Zfc/0.01 Zq) . (25) 

A similar calculation can be made for the first two rota- 
tional lines of H2, at 28.22 fim and 17.03 /im, assuming 
LTE populations for the upper rotational levels, J = 2 
and 3. This calculation differs somewhat from the numer- 
ical calculations, in which we fix the ortho/para ratio at 
a constant value 3:1. For the 2-level analytic model, 

L28t.m = (1.4 X 10^7 erg s-^)Ms (/H2/O.OOI) (26) 
iiT;.™ = (1.2 X 10^8 erg s-^)Ms (/H2/O.OOI) (27) 

Figure lTUI shows the individual line luminosities from a 
numerical model with 10* Mq of cooling gas at 10^^ Zq 
and residual molecular fraction /h2 = 2x 10~^ at temper- 
atures ranging from 10-1000 K. Note that the metal fine- 
structure lines of [O I] , [Si II] , and [Fe II] are considerably 



11 




o 



iogio(VK) 

Fig. 10. — Individual line luminosities for 10* Mq of gas at var- 
ious temperatures, for fine-structure lines of [C II] (158 /^m), [O I] 
(63 fim), [Si II] (34.8 /zm), [Fe II] (25.99 /zm, 35.35 urn), and the 
lowest rotational lines of H2 at 28.22 /im and 17.03 fim. These 
models assume maximally efficient cooling, with high-density LTE 
populations, a single metallicity Z = 0.01 Zq, and molecular frac- 
tion /h2 = 2 X 10"*. 



stronger than the H2 rotational hnes. Even for metaUic- 
ities {Zo « 0.001) and molecular fraction (/h2 ~ 0.001), 
the [O I] fine-structure line is 200 times stronger than 
the 17.03 /im line of H2 and 1000 times stronger than 
the 28.22 fim line. 

At redshifts z Ki 4, the source luminosity distances are 
dL « {c/Ho)[z{l + z/4)] « (10^^ cm)h^Q. Luminosities 
of (10"*^ erg s^^)L4i therefore correspond to line fluxes 
(8 X 10~^^ W m~^) L41, where we have converted from 
cgs units to SI units of line sensitivity (W m~^) used by 
the SAFIR studies. At z « 10, is 4 times larger, so 
these fluxes would drop by another factor of 16. 

For a given metallicity, Z, the strongest lines from 
metal-enriched clouds are probably [O I] 63.18 /im and 
[Fe II] 25.99 fim. Lines of [Si II] (34.8 /im) and [C II] 
(158 /tm) are ten times weaker, at the same metallicity. 
However, as we have noted above, it is possible that O, Si, 
and Fe metallicities could be enhanced in regions contam- 
inated by products of massive stars. For 10^ Mq of cool- 
ing gas at 1% solar metallicity and T « 200 K, these lines 
should have luminosities ^ 2 x 10^^ erg s"^ and fluxes 
of (1 - 2) X 10-21 W m-2 at z = 4 - 5. Ground-based 
searches for these lines in the 350 fim band would probe 
redshift intervals zqi = 4.3 — 5.0 and zpeii = 11.8 — 13.6. 
Searches at other redshifts open up in the space FIR. 



4. DISCUSSION AND CONCLUSIONS 

In this paper, we have examined the thermal evolu- 
tion of a primordial gas that has been enriched with 
trace amounts of metals. We have assumed that 
these metals have managed to leave their parent halos 
ifrumlinson ct a l. 2004), where the first stars have ex- 
ploded as supernovae. Therefore, a maximum level of 



metallicity is reached in the IGIVI. The purpose of this 
paper was to calculate the fraction of metals necessary 
to produce continuum fragmentation and to examine the 
properties of the gas when this occurs. 

The results of this paper have been divided in two 
parts, in the first part (5 13.1(1 . we assumed a single metal- 
licity, with C, O, Si, and Fe appearing in relative solar 
abundances. We used this metallicity as a exploratory 
parameter to find the instant at which the gas frag- 
ments. The second part of our study considered indi- 
vidual metallicities (S 13. 2() . This is a more natural ap- 
proach, in which the fractional abundance for each metal 
comes from theoretical PISN y ields, averaged over var i- 
ous mass ranges, calculated bv lTumlinson et" 
In this case, we used Fe as a exploratory parameter and 
the amount of metallicity for each of the other three met- 
als considered in this paper depends on the mass range 
of the PISN. 

We have neglected in this work any kind of radia- 
tive feedback (UV or X-ray), more specifically photo- 
ionization heating, which raises the temperature of the 
gas and therefore shutting off star-formation (the tem- 
perature raises over > 10'* K). Another neglected heat- 
ing process is cosmic rays, which will be examined in a 
later paper. To compare our results with work done by 
previous authors, we looked for the point in clump evo- 
lution when the cooling rate (Acii,oi,Siii,Fcii) equals the 
adiabatic co mpressional heating rat e. This is the defini- 
tion used bv lBromm fc Loebl tOO^ for [C II] and [O I] 
cooling. We have also introduced fine-structure cooling 
by [Si II] and [Fe II], and we explored the slight differ- 
ences when one equates the total cooling rate (including 
molecular cooling) to the compressional heating rate. 

Our results at high densities agree well with 

iBromm fc Loebl l|2003j) . However, the required level of 
metallicity for a gas to fragment at high redshift depends 
on whether the second-generation star formation is trig- 
gered from clumps or from low density gas within the 
halo. As can be seen from Figures 4 and 8, the critical 
metallicity is considerably higher 1% solar) for the 
mean gas density, nu 0.30 cm"'^, typical of 10^ Mq 
virialized halos at z « 20. These halos may have denser 
clumps, into which the first heavy elements are incorpo- 
rated, but the details of these processes are complicated 
(Tumlinson, Venkatesan, & ShuU 2004). Numerical sim- 
ulations are probably required to fully understand the 
metal transport and fraction of metals that escape into 
the lower-density IGM. 

At densities n > Tier for the important fine-structure 
coolants, the fragmentation criterion shifts ZcHt to higher 
values. Cooling is most efficient at or near the densities, 
n « ricr, at which the fine-structure levels reach LTE. For 
this density environment, the minimum critical metallic- 
ities, Zerit, are [C/H]erit ~ -3.48, [0/H]erit « -3.78, 

[Si/H]crit = -3.54, and [Fe/H]crit ~ -3.52. When 
all metal transitions are included, we find a minimum 
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Zcrit ^10^'^-°^ Zq and Afj « 117 A/q, with a tempera- 
ture Tciit = 199 K. In this case, aU the metals appear in 
the cooHng function with the same fractional abundance, 
but these results are highly dependent on the amount 
of mixing in the gas, which may determine the relative 
abundances of each metal. 

An important parameter to consider in more detail 
is the initial abundance of residual molecular hydro- 
gen. In all the simulations presented in this paper, the 
initial fractional abundance was taken as the residual 
II2 fraction, /H2aos = n-H2,ros/n-tot = 2.0 x 10^"*. We 
also ran th e single-m ctallicity simulations with /H2,ros = 
1.1 X 10"^ l|Cxalli fc Pallaa9981 The differences between 
simulations run with I0W-H2 and high-H2 residual abun- 
dances are small at low densities. The maximum differ- 
ence in .Zcrit is '~ 50% at ri/ « 2000 cm"-^. 

For excitation from the ground vibrational state of 
H2, the most important rotational cooling lines are at 
28.22 ^m (J = 2 ^ 0) and 17.03 urn {J = 3 ^ 1), 
labeled (0-0) S(0) and (0-0) S(l), respectively. Because 
of their low radiative transition rates, these II2 emission 
lines are quite weak; the metal fine-structure lines are ex- 
pected to be considerably stronger, even at T = 200 K, 
once the metallicity rises above 10"'* Zq. 

Using individual metallicities for C, O, Si, and Fe, 



we found an important distinction from the single- 
metallicity case. These results suggest that, for massive 
stars enhanced in a-process elements (O, Si and Fe), the 
critical metallicity for fragmentation is much higher than 
previous works have suggested. Higher metallicities are 
also required if the metal-enriched gas has lower densities 
than rtcr- 

Redshifted fine-structure lines from metal-enriched pri- 
mordial gas at 200 K could be observable from high- 
redshift halos. In addition to the weak rotational lines 
of H2 (28 Aim and 17 ^im) and [C II] 157.74 ^m, the 
most important coolants are likely to be [O I] 63.18 /im, 
[Si II] 34.8 ^m, and [Fe II] (25.99 ^lm and 35.35 nm). 
With redshifting, these lines would appear in the FIR 
sub-mm, and even the mm-band, if (1 + ^) w 10 — 20. 
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TABLE A5 
Chemical Reactions and Rates^ 



Reaction Rate (cm"^ s ^) Reference 



1 


H + 


e — 


■> H+ + 2e- 


5.9 X 10-"T(5''5(l + r|''5)-iexp(-1.58/T5) 


HTL 


2 


H+ - 


f e- - 




3.3 X 10-10 T(,-'^'^(l + r0'7)-i 


HTL 


Q 
O 


XT "h 


e — 


-> n -f- 7 


1 zl V in— 18 'pO.QS ovi-if' 7^, /I f^9^ 
1.4 X lU Q exp^^ — _Z4y'±.DZ^ 




A 

4 


TJ — 
t± 


^ tl — 


H2 + e 


i.o X iU 




5 






> Zri 


/I n \y in — 6 — 0.5 
4.0 X 10 J 


r O 


6 


H + 


H+ - 


-H++7 


2.1 X 10-23 T(,l 8exp(-2/Ti) 


SLD 


7 


■r-r-l- 


f H — 


Hz + H+ 


6.4 X 10-1" 


GP 


8 




f e- - 


2H 


1.2 X lO-'^ Tp-" '' 


SLD 


9 


Ha H 


~H+ - 


-.H++H 


min 1^3.0 X IQ-i" cxp(-2.11/T4), 1.5 X lO-l" 
cxp(-1.40/r4)] 


GP 


10 


Hz H 


-H — 


-> 3H 


7.1 X 10-i''t2'01(i + 2.13T5)-3-'''1 
exp(-5.18/T4) 


AAZN 


11 


Hz H 


- e~ — 


2H + e- 


4.4 X 10-"'T(5''35exp(-1.02/T5) 


GP 


12 


He 4 


- — 


He+ + 2e- 


2.4 X 10-11 (1 + T0'5)-i cxp(-2.85/r5) 


HTL 


13 


He+ 


+ e- 


> He + 7 


1.5 X 10-1" T-°-^^ + 1.9 X 10-3 Tg-i '^ 
exp{-5.64/T5) (0.3 + exp(9.40/T4)) 


AAZN 


14 


He+ 


+ e- 


> He+2 _^ 


5.7 X 10-12 T°-^ (1 + TgO-S)-! cxp{-6.32/T5) 


HTL 


15 


He+2 + e- 


> He+ + 7 


1.3 X 10-9T-"'^(1 + T0'7)-i 


HTL 



''This table summarizes the important reactions needed to calculate accurately the abundances of Hz and e . References are: HTL 
Haiman, Thoul & Loeb (1996); GP GaUi & Palla (1998); FC Fuller & Couchman (2000); SLD Stancil, Lcpp & Dalgarno (1998); AAZN 
Abel, Anninos, Zhang & Norman (1997). 



APPENDIX 

PRIMORDIAL CHEMICAL REACTION NETWORK 

Table IA5I shows the main chemical reaction between elements in the primordial gas used in our simulations. The 
gas is co mposed by H2, H, H"* ", Hij", H~, He, He+, He^^, and e~. The rea. c tion r a tes were comp i led fr om t he following 
autho rs: iHaiman et all l|l996|) . iGaUi fc Pallal l|l998() . iFuller fc Couchma n (20001. Sta ncil et alJ l)l998(l and lAbel et"all 
l)l997|) . The minimal model shown in the table was defined bv lHutchines et aL ((2002) ■ 



